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ABSTRACT 

In order to understand the climate on terrestrial planets orbiting nearby Sun-like stars, one would 
like to know their thermal inertia. We use a global climate model to simulate the thermal phase vari- 
ations of Earth-analogs and test whether these data could distinguish between planets with different 
heat storage and heat transport characteristics. In particular, we consider a temperate climate with 
polar ice caps (like modern Earth) , and a snowball state where the oceans are globally covered in ice. 
We first quantitatively study the periodic radiative forcing from, and climatic response to, rotation, 
obliquity, and eccentricity. Orbital eccentricity and seasonal changes in albedo cause variations in the 
global-mean absorbed flux. The responses of the two climates to these global seasons indicate that the 
temperate planet has 3 x the bulk heat capacity of the snowball planet due to the presence of liquid 
water oceans. The temperate obliquity seasons are weaker than one would expect based on thermal 
inertia alone; this is due to cross-equatorial oceanic and atmospheric energy transport. Thermal iner- 
tia and cross-equatorial heat transport have qualitatively different effects on obliquity seasons, insofar 
as heat transport tends to reduce seasonal amplitude without inducing a phase lag. For an Earth-like 
planet, however, this effect is masked by the mixing of signals from low thermal inertia regions (sea 
ice and land) with that from high thermal inertia regions (oceans), which also produces a damped 
response with small phase lag. We then simulate thermal lightcurves as they would appear to a high- 
contrast imaging mission (TPF-I/Darwin). In order of importance to the present simulations, which 
use modern-Earth orbital parameters, the three drivers of thermal phase variations are 1) obliquity 
seasons, 2) diurnal cycle, and 3) global seasons. Obliquity seasons are the dominant source of phase 
variations for most viewing angles. A pole-on observer would measure peak-to-trough amplitudes of 
13% and 47% for the temperate and snowball climates, respectively. Diurnal heating is important 
for equatorial observers (~ 5% phase variations), because the obliquity effects cancel to first order 
from that vantage. Finally, we compare the prospects of optical vs. thermal direct imaging missions 
for constraining the climate on exoplanets and conclude that while zero- and one-dimensional models 
are best served by thermal measurements, second-order models accounting for seasons and planetary 
thermal inertia would require both optical and thermal observations. 



1. INTRODUCTION 

An extrasolar planet's disk-averaged emitting temper- 
ature changes throughout its orbit, a phenomenon known 
as thermal phase variations. In this paper we use a com- 
prehensive global climate model (GCM) to compare the 
thermal phase variations for two very different climates 
that are known to have existed in the Sun's habitable 
zone (HZ): a temperate climate, and a pan-glacial cli- 
mate from Earth's geological past ("snowball Earth"). 
After introducing our climate model (Section 2) and dis- 
cussing its time-averaged characteristics (Section 3), we 
study the periodic radiative forcings and intrinsic climate 
response (Section 4). We then consider, in Section 5, how 
these climate variations manifest themselves to an exter- 
nal observer; we discuss our results in Section 6. Such a 
study can help determine whether time-resolved thermal 
photometry will be as useful a tool for directly imaged 
exoplanets as it is for the current crop of short-period 
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giant planets. 

As iCowan et al.l (|2011| ) noted, the most direct way to 
distinguish a temperate climate from a snowball is to 
measure the planet's Bond albedo. By combining an ex- 
oplanet's emitting temperature with its bolometric flux, 
one can obtain estimates of its radius and albedo. The 
albedo of a snowball planet would be more than twice 
that of a terrestrial planet, making such a measure- 
ment easy for an infrared direct imaging mission (e.g., 
NASA's Terrestri al Planet Finder Interferometer TPF-I, 
or ESA's Darwin: iBeichmanll 19981: lFridlundll2000l respec- 
tively) However, an albedo estimate is necessary, but 
insufficient, for understanding a planet's climate. Per- 
versely, the end state at the inner edge of the HZ is also 
highly reflective, if Venus is any indication. 

A climate constraint orthogonal to albedo is plane- 
tary thermal inertia, which depends on atmospheric and 
surface characteristics. Terrestrial planets may be di- 
vided into four categories based on albedo and thermal 
inertia: 1) low albedo and low thermal inertia (Mars), 

5 It is worth noting, however, that the albedo of snow, ice 
and clouds is much less for planets orbiting low-mass stars 
llJoshi fc Habcrlc 2012). This makes it harder for such planets 
to enter a runaway snowball, but also makes it harder to recognize 
such snowballs, once formed. 
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2) high albedo and low thermal inertia (snowball), 3) 
high albedo and high thermal inertia (Venus), and 4) 
low albedo and high thermal inertia (Earth). Further- 
more, a large thermal inertia for a planet with a modest 
atmosphere would constitute circumstantial evidence for 
liquid water oceans. 

A planet's thermal inertia can be estimated 
based on its resp o nse to external radiative forcing. 
iGaidos fc Williamsl f|2004T ) used an analytic energy 
balance model (EBM), which assumes diffusive heat 
transport, to explore how the thermal phase variations 
of terrestrial planets are affected by a planet's eccentric- 
ity, obliquity and thermal inertia. They found strong 
degeneracies between these factors: even if the planet's 
orbit is well characterized via astrometry or imaging, 
the obliquity, equinox phase and thermal inertia cannot 
be uniquely determined without detecting second-order 
phase variations. 

The present study goes beyond the important work of 
IGaidos fc W illiams (l200l in the following ways: 1) we 
use a GCM rather than an analytic EBM, allowing us 
to self-consistently model the three-dimensional atmo- 
sphere and ocean dynamics of these planets, 2) we test 
the effect of planetary thermal inertia using temperate 
and snowball climates, which are realistic and stable on 
geological timescales, 3) we include the effects of diur- 
nal heating and rotational variability, and 4) we simu- 
late ten years worth of observations for each planet, al- 
lowing us to evaluate the repeatability of these patterns 
given stochastic weather and inter-a nnual variability. Fi- 
nally, our study is complementary to IGaidos fc Williamsl 
(2004): we are unable to vary eccentricity or obliquity be- 
cause of computational expense, but we consider a wide 
variety of viewing geometries. 

1.1. Viewing Geometry 

It is useful to first define the system geometry, shown 
in Figure [TJ An exoplanet's orbital inclination, i, denotes 
how far it is from being in a face-on orbit: a planet with 
i = orbits in the plane of the sky, while i = 90° means 
an edge-on orbit. Unless a planet is in a perfectly face- 
on orbit, it exhibits phases: when the planet is closest 
to us, inferior conjunction, we primarily see the planet's 
nightside; when the planet is farthest from us, superior 
conjunction, we see mostly dayside. A planet on an edge- 
on orbit is "new" at inferior conjunction, and "full" at 
superior conjunction. Quadrature occurs when we see 
half of the illuminated hemisphere of the planet (as with 
the Moon's quarter phases). A planet on a face-on or- 
bit is permanently in quadrature; for other inclinations, 
quadrature occurs twice per orbit. 

The angle between the planet's rotation axis, LJ ro t, and 
its orbital axis, w rb is the obliquity, 0. The orientation 
of the rotational axis with respect to an observer is spec- 
ified by the equinox phase, £e- Note that equinox phase 
is related to — but not the same as — precession phase, 
the orientation of the planet's spin axis with respect to 
periastron. 

The sub-observer latitude, 6> b s , is the latitude of the 
point directly beneath the observer. If the planet's equa- 
tor is directly beneath the observer, then obs = 0° 
( "equator-on" viewing geometry), while o bs = 90° in- 
dicates that the observer is located directly above the 
planet's north pole ("pole-on" viewing geometry). 




Plane ' 

Fig. 1. — The viewing geometry is defined by any three of i, 
dobs, Q, or £e- Orbital inclination, i is the angle between the 
orbital plane of the planet and the sky plane, or equivalently 
the angle between the orbital angular momentum and the 
line of sight. The sub-observer latitude, is the latitude 

directly below the observer. These two angles are related via 
the obliquity, 0, and equinox phase, £e, via Equation [T] 

The sub-observer latitude is related to the inclination, 
planetary obliquity, and equinox phase: 

sin (9 obs = sin cos £e sin i + cos cos i, (1) 

where conservation of angular momentum dictates that 
all of the parameters in the equation are constants, pro- 
vided that we neglect precession. 

1.2. Sources of Thermal Phase Variations 

An externally heated planetary atmosphere can 
roughly be thought of as a heat engine. For example, 
surface-to-tropopause and equator-to-pole heating gra- 
dients on Earth lead to steady-state transport of energy: 
upward and poleward, respectively. It is the higher-order 
periodic climate forcings, however, that cause thermal 
phase variations (Table Q] and Figure [2]). These varia- 
tions are damped by heat storage (i.e., thermal inertia) 
and transport (atmospheric and/or oceanic). 

Diurnal Cycle: A planet receives flux from its host 
star during the day but not at night. This tends to make 
a planet's day-side warmer than its night side, leading 
to phase variations: the day-side is most visible near 
superior conjunction, and the night side is most visible 
near inferior conjunction. Insofar as a planet has some 
thermal inertia, the day-night temperature contrast is 
reduced and the hottest time of day occurs in the af- 
ternoon rather than at noon. If the advective timescale 
is comparable to, or shorter than, the rotational period, 
then zonal (E-W) heat transport becomes an important 
factor for damping the diurnal temperature swings. The 
diurnal cycle is most observable for planets on edge-on 
orbits, and is the sole source of phase variations for a 
zero-eccentricity, zero-obliquity planet. 

If a planet's surface has inhomogeneous emitting tem- 
perature (e.g., due to patchy clouds), then different re- 
gions of the planet will have different emitting tempera- 
tures at the same local solar time, leading to rotational 
variability. This manifests itself to an external, inertial 
observer as high-frequency variations in emitting tem- 
perature, where the magnitude of variations depends ul- 
timately on the longitudinal albedo and/or thermal iner- 
tia contrasts on the planet. While this effect is included 
in our simulations and could in principle be used to deter - 
mine a planet's rotation rate ([Gomez-Leal et aTI ?2012'). 
it would be averaged over for realistic integration times 
(> 24 hrs), and we do not discuss it further in this paper. 
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Sources of Thermal Phases 



Apastron: 

Minimum Incident Flux 



Superior Conjunction: 
Day-Side Visible 




nferior Conjunction: 
Night-Side Visible 




Winter Hemisphere Visible 



Fig. 2. — The three periodic drivers of climate, the diurnal cycle, obliquity seasons, and global seasons, are the sources of 
thermal phase variations. The response of a planet to these drivers depends on its thermal inertia and heat transport, as 
described in the text and in Table [l] 



TABLE 1 

External Periodic Climate Forcing and the Sources of Thermal Phases 



Climate Forcing 


Frequency 


Amplitude" 1 


Damping Mechanism 


Favored Geometry 


Diurnal Cycle 
Eccentricity Seasons 
Obliquity Seasons 


Rotational 6 

Orbital 

Orbital 


oc S (l - A) 

„. 4eS (l-A) 

oc S (l - A) sinG 


Heat Storage and Zonal Transport 11 
Heat Storage 

Heat Storage and Meridional Transport 


Edge-On Orbit 
Always Present 
Pole-On Rotation 



a The intrinsic pcak-to-trough forcing amplitude, S max ~~ 5'min. is expressed in terms of the Solar constant, Sq, the planet's 
Bond albedo, A, orbital eccentricity, e, and obliquity, 0. The more familiar relative amplitude for eccentricity seasons is 
Smax/Smm = (f 1 + e ) / ' i 1 ~ e )) 2 • Thc sin scaling for obliquity seasons is approximately correct for mid— high latitudes and 
thus for thc hemispherical average (we determined this numerically solving thc expressions in Section 7.3 of Picrrchumbcrt 2010). 

The diurnal forcing occurs at the rotational frequency, but it manifests itself to a distant observer on the orbital timcscale as 
phase variations. 

c On Earth, wind speeds arc much slower than rotational velocity for all but polar regions, so thc diurnal cycle is almost entirely 
damped by thermal inertia. 



Eccentricity Seasons: A planet with non-zero orbital 
eccentricity tends to be hottest near periastron and 
coolest near apastron. If the planet has some thermal in- 
ertia, the seasonal temperature swings are damped, and 
the maximum temperature occurs somewhat after pe- 
riastron. Orbital eccentricity affects the entire planet's 
power budget, and hence affects thermal phase variations 
seen from any vantage point; eccentricity effects may be 
isolated if a zero-obliquity planet is viewed face-on. 

Obliquity Seasons: Non-zero obliquity leads to the 
northern hemisphere being hotter during northern sum- 
mer and cooler during northern winter. Obliquity sea- 
sons are damped by thermal inertia and meridional (N-S) 
heat transport. Insofar as one hemisphere is more vis- 
ible to the observer (# obs ^ 0), obliquity seasons cause 
thermal phase variations. If the northern and southern 
hemispheres have very different albedos and/or thermal 
inertia, the planet's global power budget will exhibit sea- 
sonal variations regardless of viewing geometry, much as 
with eccentricity seasons. In general, obliquity seasons 
are most extreme as seen by a pole-on observer. A planet 
on a circular, face-on orbit will have thermal variations 
entirely dictated by obliquity seasons. 

1.3. Snowball Climate 



A t errestrial planet in the HZ of its host star (|Hartl 
ITffllAb^[l99aiKasting et al.lll993ft does not necessar- 
ily have a temperate climate or harbor liquid water on 
its surface. The high reflectivity of ice and snow com- 
pared to land and ocean allows an Earth-like planet at 
the exact same distance from its star and with the same 
atmosphere to exist in either a temp erate state, like mod- 
ern Earth, or an icy snowball st ate (jBudvko 1969; Sellersl 
[19691 IMarotzke fc Botzetll20fll . 

It is believed, for example, that Earth has on occasion 
(~710 Ma, ^635 Ma, and possibly at other times) fallen 
into a s nowball state, during which global oc eans com- 
pletely dKirschvink|]1992t Hoff man et all 1 1998ft or nearly 
completely (|Hvde et al.ri200Ct lAbbot et all 120 lift froze 
over. As a limiting case, we will be considering here 
a planet with entirely ice-covered oceans. The proxi- 
mate mechanism for the initiation of a snowball state 
is a runaway ice/snow- albedo feedback: replacing ocean 
with sea ice leads to cooli ng and more ice formation 
()Budvkolll969l: ISellersI 1 1969ft . The ultimate cause of the 
specific snowball events in Earth 's history is still un- 
known (|Pierrehumbert et al.ll201lft . 

The high albedo of ice and snow make a snowba l l state, 
once formed, extremely stable (jBudykd Il969t ISellersI 
Il969t IMarotzke fc Botzetl 120071 ). Compared to Earth's 
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4.5 Gyr geological history, snowball events have been 
rare and brief (millions of years), with a duty cycle of 
order 1%, probably du e to the geologically rapid silicate- 
weathering feedback (jWalker et alJ Il98lh . The glob- 
ally frozen continents and low temperatures halted sil- 
icate weathering and the associated pulldown of atmo- 
spheric CO2, but greenhouse gases kept being resup- 
plied by volcanism and plate tectonics. The resulting 
increase in greenhouse warming — in conjunction with 
other effects such as the dirt ying of tropical sea ice 
(jAbbot fc Pierrehum bert 201(J) — snapped Earth out of 
the snowball state. 

It is not known how the escape from a snowball state 
would play out on exoplanets. Possible differences in- 
clude the presence of a large water reservoir (wet plan- 
ets are more prone to r unaway snowball than "Dune" 
planets; lAbe et ail 1201 incident stellar flux (reduced 
power budgets make planets more susceptible to global 
glaciations, including th e limit of a starless planet, Taiika 
120081: lAbbot fc Switzerl 120111) , host star spectral type 
(bluer stars lead to a stronge r temperature-albedo feed- 
back; IJoshi fc Haberld I2012D , planetary obliquity (low- 
obliquit y planets are slightly mor e prone to snowball 
events: IWilliams fc Kastind 119971: IWilliams fc Pollard! 
120031 iSpiegel et al.ll2009D . orbital eccentricity (larger ec- 
centricity slightly increases t he propensity of a planet 
to enter a snowba ll state: IWilliams fc Pollardl 120021 : 
IDressing et al.|[2~010f) , planetary rotation and size (slower 
rotation or a smaller planet increases meridional trans- 
port, making it easier to enter a snowball state), and the 
effectiveness of an extrasolar carbon-silicate cycle (con- 
tinent configuration, water levels, nature of mantle con- 
vection, etc.). 

Most importantly for our purposes, a snowball climate 
has relatively little surface thermal inertia compared to a 
temperate climate: a few meters equivalent depth of wa- 
ter compared to the current temperate climate with its 
tens of meters equivalent depth of water Trying to dis- 
tinguish between these climates based solely on thermal 
phase variations is therefore a good test for the utility of 
thermal phase variations in constraining the heat capac- 
ity of HZ terrestrial exoplanets. 

2. GLOBAL CLIMATE MODELS 

To simulate thermal phase variations for a variety of 
viewing geometries, we use output from the coupled 
ocean-atmosphere global climate model ECHAM5 /MPI- 
OM. A detailed description of the GCM is given in 
IVoigt et all f|201lD and references therein. We use top-of- 
atmosphere (TOA) incident stellar flux, outgoing long- 
wave radiation, and reflected light, each of which is out- 
putted every two hours at each grid point. 

Specifical l y, we here use two simulations presented in 
IVoigt et al.l (| 2 1 If ) . These simulations use a modern- 
Earth solar constant of 1367 Wm~ 2 , pre- industrial atmo- 
spheric carbon dioxide (278 ppm), a modern-Earth orbit 
(23.5° obliquity, 1.7% eccentricity) and two large tropi- 
cal supercontinents that are meant to represent Earth's 
topography 635 Ma before present. The time-averaged 
albedo and temperature maps of the two models are 
shown in Figure [3] 

6 Earth's oceans are on average 3 km deep, but only the top tens 
of meters participate in damping seasonal variations. 



The "temperate" simulation is initiated from a climate 
state with polar ice caps and remains in such a state 
throughout. The "snowball" simulation is initiated from 
a climate state with globally ice-covered oceans and re- 
mains in such a state throughout. As found in previ- 
ous snowball climate simulations, the low latitudes, and 
in particular the majority of the continental area, re- 
main bare because of a reduced hydrologi cal cycle and 
net evaporation (e.g.. iPierrehum bcrt 2005). Both simu- 
lations were run to equilibrium, then for an additional 
ten years, with shortwave and longwave radiation out- 
puts saved every two hours. 

Although our temperate simulation is roughly com- 
mensurate with the modern-day Earth, it differs in two 
important ways: 1) the continental configuration is much 
more N-S symmetric, reducing the impact of obliquity 
on the global power budget, and 2) the continents are 
more uniform longitudi nally, reducing the ampl itude of 
rotational variations (cf. iGomez- Leal et al.|[20ia) . These 
differences make the resulting climate variations of our 
temperate planet, if anything, more straightforward to 
interpret than the modern Earth, strengthening our ul- 
timate conclusions. 

3. TIME-AVERAGED POWER BUDGET 

The time-averaged, total absorbed power of the planet 
must equal its emitted thermal power. We verify that 
this is indeed the case, to one part in a thousand, for the 
temperate simulation, while the snowball's emitted flux 
exceeds its absorbed flux by 5.5%. More energy leaves 
the atmosphere at the top than enters at the bottom, 
yet the snowball planet does not cool off, implying that 
there is a small spurious heat source in the model atmo- 
sphere. This discrepancy, which may be a surface diag- 
nostic bug rather than a true model error, was noted by 
IVoigt et al.l (|2011| ) and should not significantly affect our 
results, since we only compare time variations in ther- 
mal emission. We simply renormalize the snowball simu- 
lation's thermal emission to the time-averaged absorbed 
power in Sections 3 & 4 of this paper. 

3.1. Year-Averaged Meridional Heat Transport 

The modest obliquity of our simulated planets dic- 
tates that high latitudes receive less orbit-averaged in- 
solation than low latitudes. In Figure [4] we show the 
top-of- atmosphere radiative budget for the two simula- 
tions. Averaging over the seasonal variations that will 
be the subject of the bulk of this paper, the equator- 
pole gradient in absorbed flux (black lines) results in a 
meridional temperature gradient, which drives poleward 
energy transport. As a result, the equator-pole gradient 
in thermal emission (red lines) is muted compared to the 
gradient in absorbed flux. 

Comparing the radiative budgets of the two simula- 
tions, it is clear that the meridional radiative forcing 
gradient is somewhat greater for the temperate case than 
for the snowball (240 vs. 200 W/m 2 ), while the resulting 
gradient in thermal emission is similar for the two planets 
(80 vs. 70 W/m 2 ). This indicates that more poleward 
heat transport is occurring on the temperate planet. 

Figure [5] shows the orbit-averaged northward heat 
transport for the two simulations; heat transport on the 
temperate planet is roughly twice that on the snowball. 
Some of this difference can be attributed to the weaker 
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Fig. 3. — Ten-year averaged maps of top-of-atmosphere Bond albedo (top panels) and effective temperature (bottom panels) for the 
temperate (left) and snowball (right) models. Green lines indicate the coasts of the equatorial continents; the blue lines show the annual 
mean extent of sea ice in the temperate model. 



meridional forcing on the snowball planet. Moreover, 
transport is much greater in the temperate climate be- 
cause of two compounding effects: the presence of ocean 
transport and greater atmospheric humidity. 

On their own, neither the absence of ocean heat trans- 
port, nor the drier air in principle need to grossly af- 
fect the total meridional heat transport of an Earth- 
like planet due t o "compensations" noted an d ex- 
plained bv IStonel (119781) . IVallis fc Farnetil (pnM ). and 
lEndertonfc Marshall (2009?). But in combination, these 
two changes make it impossible for the climate to fully 
compensate, with the net effect that the mid-latitude 
meridional heat transport in the snowball simulation is 
significantly reduced. 

4. PERIODIC CLIMATE FORCING 

We now consider periodic changes in the planet's power 
budget. We first develop a simple model for interpreting 
time- variations in absorbed and emitted power, then con- 
sider each of the periodic climate forcings in turn, as they 
manifest themselves in our simulations. The planet's re- 
sponse to radiative forcing is independent of an external 
observer, but is the ultimate cause of the thermal phase 
variations discussed in Section 5. 

4.1. Response to Periodic Forcing 



In order to quantitatively estimate thermal inertia 
from periodic variations in absorbed and emitted power, 
we approximate the clim ate system with the following 
differential equation (e.g.. iPierrehum bcrt 20l3): 



dT 



A^abs COs(wi) - bT, 



(2) 



where /i is the effective thermal inertia coefficient, T is 
the emitting temperature, Ai^bs is the semi-amplitude 
of variations in the absorbed power, lo is the driving 
frequency, and b linearly approximates the relation be- 
tween temperature and emitted flux (the exact differen- 
tial equation, analytic a nd numerical solution s are devel- 
oped for a blackbody in lCowan fc Agolll2011l ). 

It is worth noting that this formalism is usually 
adopted for surface temperature variations. Nonethe- 
less, the most accessible temperature for any planet other 
than Earth is the emitting temperature so it is natural 
to consider effective thermal inertia as it affects top-of- 
atmosphere absorbed and emitted flux. Inferring bulk 
thermal inertia for the underlying atmosphere and sur- 
face is of course complicated when one is only sensitive 
to the emitting temperature. For example, the role of 
greenhouse gases in Equation [5] is normally to reduce 
b, since for a given change in surface temperature, less 
flux escapes to space. In the current scenario where 
T denotes top-of-atmosphere emitting temperature, the 




Fig. 4. — Top-of-atmosphcre absorbed (black) and emitted (red) flux as a function of latitude. 




lotitude 



Fig. 5. — Seasonally-averaged northward heat transport for the 
temperate (black) and snowball (blue) climates. The snowball 
planet has ice-covered oceans and thus negligible oceanic heat 
transport. 

emitted flux is erT 4 , where a is the Stefan-Boltzmann 
constant. In the limit of small temperature excursions, 
we may adopt b = 4<tTq , where T is the average emit- 
ting temperature. The greenhouse effect is therefore 
folded into /i: greater infrared opacity increases the at- 
mosphere's radiative timescale and hence the lag time be- 
tween radiative forcing at the surface and response at the 
top-of- atmosphere. Insofar as upward energy transport 
in the troposphere is primarily accomplished via convec- 
tion, the atmosphere's infrared opacity determines the 
pressure of the emitting layer. Either way, the net result 
of increasing the concentration of infrared absorbers is 
increasing the effective thermal inertia coefficient. 
The periodic solution to Equation [2] is 

T(t) = ^+ ATcos(wt-$), (3) 
o 

where the temperature oscillations have a semi- 
amplitude AT = AFabs/ (byl + u) 2 t 2 ), the phase offset 
is given by tan <f> = wr, and the thermal relaxation time 
is t = n/b. 

In this linearized form, the emitted flux is proportional 
to temperature, AF cm ; = &AT, so one can estimate the 



radiative relaxation time based on the amplitudes of vari- 
ations in absorbed and emitted power: 

4.2. Global-Mean Temperature Variations 

The incoming shortwave radiation undergoes a 6.8% 
peak-to-trough orbital modulation because of Earth's 
1.7% eccentricity, while the albedo exhibits seasonal vari- 
ations of 6% and 4% in the temperate and snowball cli- 
mates, respectively. The albedo fluctuations are primar- 
ily driven by meridional (N-S) differences in albedo cou- 
pled with seasonal variations in the sub-solar latitude 
(|Cowan et al.l 12012). Seasonal changes in cloud cover 
may also contribute to the changing albedo. This effect 
cannot be predicted a priori but is an emergent phe- 
nomenon of the GCM. 

The combination of changing incident flux at the top 
of the atmosphere and changing albedo leads to peak- 
trough variations in absorbed power of 6% and 8% for 
the temperate and snowball climates, shown with black 
lines in Figure [6] (in physical units, the two simulations 
experience approximately the same 7 PW forcing) Q The 
changes in absorbed flux drive "global" seasons. The 
critical difference between global seasons and the usual 
seasons we experience on Earth are that the seasons are 
in phase in the northern and southern hemispheres, so 
cross-equatorial heat transport cannot damp them. In- 
sofar as the typical eccentricity of exoplan ets is 0.1-0.3 
(|Ford fc Rasidl2008tlMoorhead et al.ll201lD . it is possible 
that eccentricity-driven global seasons will be the norm 
on other planets. 

The thermal inertia of the planet damps the ampli- 
tude of the oscillations and produces a phase lag, as seen 
in the planet-integrated emitted power, shown in red in 
Figure |6] (peak-trough seasonal changes of 1% and 3% 
for the temperate and snowball models, respectively). 
Note that the low-amplitude response of the temperate 
planet makes the phase offset difficult to quantify. For 
the temperate planet, peak emission occurs nearly simul- 

7 The peak-to-trough amplitudes are determined by fitting a 
sinusoid to the absorbed and emitted curves. 




Fig. 6. — Planet-integrated, top of atmosphere fluxes for a temperate climate (left) and a snowball climate (right). The plots were made 
by combining ten years of simulations and averaging over rotational variability. This plot, and most others in the paper, uses the same 
relative flux scaling (left axes) but not the same absolute scaling (right axes) for the temperate and snowball simulations. 



taneously with the peak absorption, but the minimum 
emission lags the minimum in forcing by many months. 

Using Equation|3]we obtain r = 343 days and 145 days 
for the temperate and snowball climates, respectively. 
The temperate planet therefore has a thermal relaxation 
time 2.4x longer than that for the snowball. Note that 
these are longer than the atmospheric thermal relaxation 
times because we are accounting for surface heat capac- 
ity- 
It is often convenient to express the thermal inertia co- 
efficient, fi, as an equivalent mixed layer depth (MLD) 
of water, since oceans are the dominant source of ther- 
mal inertia on Earth: H — (i/(pc p ), where p and c p are 
the density and specific heat capacity of water. Using 
the globally-averaged emitting temperatures of the two 
simulations (To = 249 K and 232 K), we estimate that 
the temperate and snowball planets have MLDs of 25 
and 8 meters, respectively; i.e., the temperate planet has 
approximately 3x the thermal inertia of the snowball 
planet. 

The greater thermal inertia of the temperate planet is 
primar ily due to the presence of surface liquid water. Fol- 
lowing [inerrehumbert (2010), Earth's oceans have an av- 
erage mixed layer depth of 20-50 m, solid surfaces (con- 
tinent or sea ice) have effective MLDs of 2-3 m, and the 
atmosphere contributes another 2-3 m. Accounting for 
the 57% coverage of surface liquid water in our temperate 
model, one expects it to have 2.4-5.5 times greater ther- 
mal inertia than the snowball planet, in good agreement 
with our simulations. 

4.3. Obliquity Seasons 

In Figure [7] we again consider total absorbed and 
emitted flux, but separate the northern and southern 
hemispheres in order to isolate the effects of obliquity. 
This approach tends to underestimate the magnitude of 
the seasons, because most of the flux is emitted from 
low-latitude regions that only experience small seasonal 
changes in incident radiation. In fact, the deep trop- 
ics receive the most insolation at equinoxes, rather than 
summer solstice. Nevertheless, the seasonal forcing and 
response are clearly visible. 

The obliquity seasons are approximately 6x greater 



for the snowball than for the temperate climate (26% vs. 
4%, as shown by the red lines in Figure [7]). If we treat the 
northern and southern hemispheres independently, ie: 
ignoring cross-equatorial heat transport, we may again 
use Equation [4] to convert the absorbed/emitted ampli- 
tudes into thermal relaxation times. We obtain values of 
r = 1044 and 155 days for the temperate and snowball 
climates, respectively. 

Since the global and obliquity seasons have the same 
forcing frequency, one expects the thermal inertia damp- 
ing for both sorts of seasons to be the same. The large 
inferred value of r for the temperate simulation obliquity 
seasons is a testament to the importance of heat trans- 
port from the summer to the winter hemispheres. Fur- 
thermore, since the obliquity seasons are so much weaker 
in the temperate simulation, one might conclude that its 
cross-equatorial heat transport is much greater than that 
of the snowball planet, but this is not the case. 

In Figure|8]we plot the seasonal cross-equatorial north- 
ward heat transport in both simulations. For the tem- 
perate planet, the ocean and atmosphere both play a 
role, while the frozen ocean in the snowball climate dic- 
tates that its heat transport is entirely atmospheric. 
Despite the lack of ocean transport, the seasonality of 
cross-equatorial energy transport is similar — in absolute 
units — for the two simulations. 

The temperate planet has a somewhat larger seasonal 
radiative forcing (40 PW vs. 30 PW; Figure!?]) and much 
greater annual-mean heat transport than the snowball 
planet (cf. Figure [5]). In the absence of thermal in- 
ertia, one would therefore expect the temperate planet 
to transport much more energy from summer to winter 
hemisphere than the snowball. But the greater thermal 
inertia of the temperate planet leads to a smaller cross- 
equatorial temperature gradient. The greater heat trans- 
port efficiency and smaller temperature gradient conspire 
to produce nearly identical cross-equatorial energy trans- 
port in the temperate and snowball simulations. 

The obliquity seasons have a much larger amplitude 
for the snowball planet than for the temperate planet, 
but the phase lag is nearly the same (1-2 weeks) in both 
cases. This is partially due to the qualitatively differ- 
ent effects of heat storage vs. heat transport in damp- 
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Fig. 7. — Absorbed (black) and emitted (red) flux for the northern (dotted) and southern (dashed) hemispheres in the temperate (left) 
and snowball (right) models. 
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Fig. 8. — Seasonal variations in the cross-equatorial heat trans- 
port. The snowball planet has globally frozen oceans and thus 
negligible oceanic heat transport. 



ing obliquity seasons: heat storage tends to increase the 
phase lag while meridional heat transport docs not (this 
is m anifestly true if heat tran sport is strictly diffusive, 
as in lGaidos fc Williams! l200l . 

The other reason that obliquity seasons can have very 
different amplitudes with similar phase lags is the inho- 
mogeneous surface on the temperate planet. This is im- 
portant for our temperate simulation because the equa- 
torial landmass and polar ice conspire to leave only 57% 
of the surface covered in liquid water. We find that 
the top-of-atmosphere thermal emission over sea ice ex- 
hibits large seasonal temperature variations that are ap- 
proximately in phase with the forcing, while that over 
ocean exhibits mu ted seasonal variations. As noted by 
IStine et all ((2009), mixing these signals results in a low 
amplitude, small-offset response, because the large am- 
plitude component (solid surfaces) dominate the offset. 
Since the land in our simulation is mostly tropical, its 
TOA thermal emission peaks during the relatively dry 
winter months, while the high surface temperature is 
masked by overlying clouds and humidity in the sum- 
mer. 



4.4. Diurnal Cycle 

To isolate the diurnal cycle, we average over the en- 
tire planet and entire orbit, considering absorbed and 
emitted flux as a function of local solar time, shown in 
Figure [9] For an Earth- like planet, zonal (E-W) heat 
transport is slow compared to planetary rotation, so the 
diurnal cycle is primarily damped by the planet's thermal 
inertia. 

We again use EquationUand obtain thermal relaxation 
times of r = 16 and 9 days for the temperate and snow- 
ball climates, respectively (or, equivalently, MLDs of 1.1 
and 0.5 meters, respectively). These heat capacities are 
much smaller than those inferred from seasonal forcing 
because much less material senses the diurnal cycle. For 
land and sea ice, downward heat transport is conductive, 
so high frequency radiative forcing does n ot penetrate as 
deep as longer-period forcing, fi cx (jPierrehumbertl 
2010). Upward heat transport in the atmosphere and 
downward heat transport in liquid water are not, strictly 
speaking, diffusive, but nonetheless only a relatively thin 
layer of the atmosphere and a shallow layer of the ocean 
respond to, and damp, the diurnal cycle. 

Furthermore, the poor vertical ice resolution in our 
GCM tends t o exaggerate the d iurnal cycle in surface 
temperatures (| Abbot et alJl2010f ) . This will most affect 
the snowball simulation, leading us to underestimate its 
thermal inertia. 

The diurnal radiative forcing peaks at noon. In the 
high-frequency forcing limit, we expect phase-lags of 7r/4 
and 7r/2 for diffusive and mixed- layer heat storage, re - 
spectively (e.g., Section 7.4.2 of iPierrehumbertl l2010f ). 
The actual phase lag is somewhere in between these val- 
ues. 



5. DISK-INTEGRATED PHOTOMETRY 

We now calculate time-resolved, disk-integrated pho- 
tometry seen by a distant observer, "light curves" . To 
compute the thermal light curves from the planet we as- 
sume diffuse emission and neglect limb-darkening: 



F=—d>F f (6,<f>,t)V(9,<t>,t)<m, (5) 
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Fig. 9. — Absorbed (black) and emitted (red) flux as a function of local solar time. The outgoing longwave radiation peaks in the late 
afternoon, while the lowest emitting temperature is reached shortly after sunrise. The angular curves are an artifact of the GCM's 2 hr 
output cadence. 

where R is the planetary radius, is the top-of- 
atmosphere outgoing longwave radiation, 9 is latitude, 
<f> is longitude, V(6, <fi, t) is the visibility of a given re- 
gion of the planet for a given observer (V is unity at the 
sub-observer point, drops as the cosine of the angle from 
the sub-observer point, and is zero on the far side of the 
planet), and § dft is the surface integral over a sphere. 

Viewing geometry only enters Equation [5] through V, 
which is a function of the sub-observer location on the 
planet. Neglecting precession, the sub-observer latitude, 
#ob s , is constant, while the sub-observer longitude evolves 
linearly in time, 4>6bs = &>bs(0) — w ro ti, where w ro t is the 
planet's rotational frequency in an inertial frame (i.e., 
the sidereal rather than solar frequency). 

For each viewing geometry we generate 10 years of sim- 
ulated thermal phase variations and track the root-mean- 
squared orbit-to-orbit variability. This inter-annual vari- 
ability is < 1%, ranging between 0.2-1.5% for most view- 
ing geometries. In Figure [10] we plot only one year of 
photometry. For consistency, we always start the obser- 
vations at inferior conjunction, when the planet is closest 
to the observer and we are seeing mostly its night-side. 

The top, middle and bottom panels of Figure [TOl show 
the same simulations, but seen by three different ob- 
servers. 

The short-term variability in the thermal lightcurve is 
due to regions with different emitting temperatures ro- 
tating in and out of view. These differences in tempera- 
ture arise because of inhomogeneities in albedo, thermal 
inertia, and diurnal cloud cycles (for a thorough discus- 
sion o f diurnal cycles in Earth's OLR, see lSmith fc Rutanl 
2003). Since we are considering integration times much 
longer than the planetary rotation rate, we do not study 
the rotational variability in this paper. 

Obliquity seasons are the largest-amplitude response 
in our simulations, so it is unsurprising that they are 
the dominant source of thermal phase variations for 
most viewing geometries. The obliquity-driven phases 
are most extreme for a pole-on observer (top panels of 
Figure I10p and most muted for an equatorial observer 
(bottom panels of Figure ITU]) . This is completely differ- 
ent from the inclination-dependance of diurnally-driven 
thermal phase variations, which are greatest for edge-on 



orbits and minimized for fa ce-on orbits (cf. iCowan et al.1 
[2lM[Crossfield et all 127)101) . 

5.1. High-Latitude or Mid-Latitude Observer 

Thermal phase variations are much greater for the 
snowball than the tem perate simulation, as one w ould ex- 
pect from Section l4~51 iGaidos fc Williama ([2004D showed 
that, in an idealized energy balance model, heat capacity 
is degenerate with obliquity in thermal phase variations. 
We cannot test this result because both of our simulated 
planets have the same obliquity, but the top four pan- 
els of Figure [TO] shows that the actual situation is likely 
worse than they concluded. The peak thermal emission 
occurs a couple of weeks after summer solstice for both 
climates, despite the muted amplitude of phase varia- 
tions in the temperate simulation, because of meridional 
heat transport and inhomogeneous surface heat capacity 
(cf. Section |4~3"]) . In other words, even if one knew the 
planetary obliquity and solstice phase independently, one 
would infer an incorrect thermal inertia from Equation[2] 

Notwithstanding the subtleties of phase lags, the am- 
plitude of thermal phase variations roughly reflect the 
intrinsic strength of obliquity seasons. The amplitudes 
of intrinsic seasons differ by a factor of 6 for the two cli- 
mates, while the obliquity-driven phase variations differ 
by 4.5. The midlatitude case (0 o h& = 45° in Figure [TO]) is 
similar to the polar observer, but with somewhat muted 
phase variations. 

The phase variations seen by a high-latitude observer 
are dominated by obliquity seasons, but the peak flux 
roughly coincides with superior conjunction in Figure [T01 
If one's intuition is based on thermal phase variations 
of hot Jupiters, it might seem that the maximum flux 
should always occur at/near superior conjunction, when 
we are seeing more of the planet's dayside. To dispel this 
notion, we plot the lightcurve for a southern observer in 
Figure [TTJ Inferior conjunction still occurs at t — 0, but 
now it is the southern seasons that affect the observed 
seasonal variability and the planet flux peaks near infe- 
rior conjunction. 

5.2. Equatorial Observer 
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Fig. 10. — Thermal photometry as seen by a pole-on (top), mid-latitude (middle), or equatorial (bottom) observer. The black vertical 
dotted line shows northern summer solstice; it is plotted for the panels with thermal phases driven by obliquity seasons. The green vertical 
dotted and dashed lines are periastron and apastron, respectively; it is plotted in the panels where the thermal phase variations are most 
affected by eccentricity. We phase-fold 9 planetary orbits worth of simulated lightcurves; the inter-annual variability is < 0.1%. The black 
points show 1-month (rather than 2hr) integrations with illustrative error bars fixed at 1%. The cartoon along the top of each panel is 
illustrative: planets that are not in edge-on orbits do not exhibit new or full phases. The viewing geometry is described by the orbital 
inclination, i (where i = is a face-on orbit and i = 90° is edge-on), sub-observer latitude, o j, a (90° at the north pole, 0° at the equator), 
and the time of periastron, t per i. 
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Fig. 11. — Obliquity-driven phase thermal variations can peak 
near inferior conjunct ion, when we are mostly seeing the planet's 
night-side. See Figure llOl for full caption. 

For the equatorial view (bottom panels of Figure [TO]) , 
the planetary obliquity should not affect the thermal 
phase variations, since 1) the northern and southern 
hemispheres are equally visible, and 2) the continents are 
nearly equatorial, so there is little to distinguish northern 
summer from southern summer. The phase variations 
from this vantage are therefore driven by global seasons 
and diurnal heating. 

If the albedo of the northern and southern hemi- 
spheres were very different, obliquity seasons would af- 
fect the phase variations seen by an equatorial ob- 
server. For example, while modern Earth's large N— 
S asymmetry in surface albedo is mostly ob scured 
by the atmosphere and clouds (iKasting et al.l 



Selsis et al.ll2007t iDonohoe fc Battistill2011b IVoigt et al.l 
20121 ). the obliquity- driven global seasons are more im- 
portant than the diurnal cycle for an equatorial observer 
i|G6mez-Leal et al.|[2012h . 

Figure [12] shows light curves for a snowball climate. 
We consider equatorial observers (6* b s = 0°) and fix the 
inclination at i = 74°, but vary the argument of peri- 
astron, which affects the amplitude and phase offset of 
thermal phases because of the co mpeting effects of glo bal 
seasons and diurnal heating (cf. iCowan fc Agol|[20rl . 

We gray-out the phases for which the planet would be 
inside the inner working angle (IWA) of a high-contrast 
imaging mission and therefore unobservable. The IWA 
depends on the details of a mission's design and is a func- 
tion of wavelength and distance. We adopt an IWA of 
34 mas and an Earth analog at 10 pc, so that phase an- 
gles within 20° of full or new phase or inaccessible (this 
is a reasonable IWA for TPF-I: lLawson et al Jl2007h . The 
planet is visible for the entire orbit if its orbital inclina- 
tion is less than 70°, as is the case for the geometries in 
Figures [TO] and HQ 

The changing power budget of the planet due to its 
non-circular orbit leads to a maximum in thermal emis- 
sion near periastron (Figure [6]). Diurnal heating, on the 
other hand, leads to a maximum in thermal emission 
between waxing quarter and superior conjunction, when 
afternoon local solar times of the planet are most visi- 
ble (Figure IH]). Note that Earth has westward equato- 
rial winds, but they are slow compared to the planet's 
rotation, so the hottest local solar time is still in the 
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Fig. 12. — Thermal phase variations seen by an equatorial ob- 
server tend to be dominated by eccentricity and diurnal effects. 
Depending on the argument of periastron (shown here as the time 
of periastron, i per j), these two effects can interfere constructively 
(top panel) or destructively (bottom panel). The grey regions de- 
note the phases that will not be accessible because the planet will 
be inside the inner working angle (IWA) of the telescope; the ge- 
ometry shown in the top panels is sufficiently c lose to face-on that 
the planet is never within the IWA. See Figure flOl for full caption. 



afternoon when averaged over the visible disk. For plan- 
ets with thicker atmospheres, the wind direction becomes 
more relevant than the sense of the planet's rotation (e.g., 
super- rotation on Venus, Titan, or hot Jupiters). 
The thermal phase variations are therefore maximized 
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when superior conjunction occurs a few months after pe- 
riastron passage, so that the diurnal and eccentricity ef- 
fects add up (top panel of Figure [12]) ; if the two effects 
are perfectly out of phase, the phase variations are at- 
tenuated (bottom panel). In all cases the phase variation 
amplitude is approximately twice as large in the snow- 
ball than the temperate simulation, as one would expect 
from Section [4] 

The diurnal phase variations are smaller than the in- 
trinsic diurnal cycle: even though a single local solar time 
is directly facing the observer at a given time, the disk- 
integrated lightcurve averages over many local times. 
Quantitatively, the first-order diurnal mode is attenuated 
by 7r/4 0.79 with respect to t he time-constant com - 
ponent for an edge-on observer (|Cowan fc Aeoll [20081. 
Since eccentricity affects the entire planet's power bud- 
get, it does not suffer from this dilution. Therefore, while 
the intrinsic diurnal cycle is nearly twice the amplitude of 
the global seasons, the diurnal contribution to the phase 
variations is only ~ 50% greater than that from eccen- 
tricity effects. 

6. DISCUSSION 
6.1. Feasibility 

Measuring thermal phase variations of an Earth ana- 
log will require space-based interferometry (e.g., NASA's 
Terrestrial Pl anet Finde r Inter f erometer or ES A's Dar- 
win mission; IBeichmanl ITqqI iFridlundl l27)00l respec- 
tively). iCockell et al.l ()2009f ) estimated that Darwin 
could obtain 20% photometry on an Earth-sized planet 
at 10 pc in 12-24 hrs of integration. Extrapolating to 
month-long integrations, one expects 2.6-3.6% photom- 
etry for an Earth-analog at 10 pc. 

We have simulated nine years of photometry, and 
found the inter-annual variability in disk-integrated flux 
is < 1%, and a mere < 0.1% after averaging over rota- 
tional variability. This opens the possibility of combining 
multiple orbits-worth of data. For example, combining 
the nine years of photometry simulated here provides a 
three-fold enhancement in the signal-to-noise ratio, or 
1% photometry with 1-month integrations. 

By using bolometric longwave radiation as a proxy for 
measured flux, we are assuming an instrument capable 
of capturing most of the planet's spectral energy distri- 
bution (SED). This is a reasonable assumption because 
one will need to measure a directly imaged planet's SED 
in order to estimate its albedo and radius. Considering 
bolometric flux also minimizes the impact of any unseen 
satellites, w hich only domin ate the system flux in narrow 
wavebands ([Robinson 20TU) . Nevertheless, the presence 
of a large satellite would considerably complicate the in- 
terpretation of thermal phase variations. 

When estimating the ability of TPF-Interferometer 
to establish a plane t's orbit and residence in the HZ, 
lLawson et al.l (|2007l ) assumed the planet has constant 
mid-infrared brightness throughout each orbit. We find 
that the thermal flux from an Earth analog exhibits or- 
bital variations of 5-50%. The amplitude of these varia- 
tions is primarily governed by the planet's obliquity and 
the sub-observer latitude, neither of which will be known 
a priori for directly imaged planets. This will compli- 
cate the retrieval of their orbital parameters from direct 
imaging. Note that the situation is yet worse at optical 



wavelengths: reflected flux variations are of order 100%, 
since the planet's night side reflects no light. 

6.2. The Utility of Thermal Phase Variations 

There is an interesting parallel between non-transiting 
hot Jupiters and directly-imaged planets (in what follows 
we assume a circular orbit, but the argument can be gen- 
eralized to non-zero eccentricity) . In the hot Jupiter sce- 
nario, the obliquity is assumed to be zero because of tidal 
effects, but orbital inclination is unknown. The thermal 
phase variations are entirely driven by diurnal heating, 
and therefore put a joint constraint on the planet's size, 
orbit al inclination, and day-night temperature contrast 
(e.g.. iCowan et al.ll2007t iCrossfield et al.ll2010[) . For di- 
rectly imaged planets, the observer latitude is unknown, 
but the orbital inclination can be measured from astrom- 
etry. For all but edge-on orbits, thermal phase variations 
can be assumed to be obliquity-driven, and therefore 
place a joint constraint on observer latitude and summer- 
winter temperature contrast. 

Focusing now on directly-imaged planets, the observer 
latitude is a function of one known quantity, the incli- 
nation, and two unknowns: the obliquity and equinox 
phase (Equation [TJ . The summer-winter temperature 
contrast is a function of obliquity, thermal inertia and 
cross-equatorial heat transport. Two measured quan- 
tities, thermal phase amplitude and offset, are clearly 
insufficient to break the four- way degeneracy^ 

Even without breaking any degeneracies, thermal 
phase variations of an imaged HZ planet might still be 
useful. Small thermal phase variations would improve 
the odds of habitability, indicating either very weak 
radiative forcing, high thermal inertia and/or vigorous 
meridional heat transport. This test does not work for 
a planet in a face-on circular orbit, which could sim- 
ply have a small obliquity. Furthermore, as stated in 
the Introduction, large thermal inertia is no guarantee of 
habitability (e.g., Venus). 

If the planet's obliquity and equinox phase were known 
through other means (see Section 16. 3p , then one could 
unambiguously estimate the summer and winter temper- 
atures of the planet. The seasonal range in temperatures 
would be a useful climate constraint in and of itself, but 
one would like to constrain the intrinsic quantities. For 
a uniform planet, one could break the two-way degen- 
eracy between thermal inertia and heat transport, since 
they have different phase offset signatures. But this logic 
does not hold for an inhomogeneous planet because the 
mixing of signals from high and low thermal inertia re- 
gions complicates the interpretation of the thermal phase 
amplitude and offset. In principle it is possible to con- 
strain inhomogeneous surfaces from rotational thermal 
variability, but this requires higher cadence and there- 
fore a larger mission. 

Insofar as the rotation rate of a terrestrial planet is 
determined by the final few major collisions of planet 
growth, it is natural to expect that some planets will 
rotate more slowly than Earth. Furthermore, planets 
on circular orbits about low-mass stars will have their 

8 The distinction between heat storage and transport may seem 
academic, but in general they are only indirectly related. On mod- 
ern Earth, for example, storage is primarily in the oceans, while 
transport is mostly accomplished by the atmosphere. 
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planetary rotation rate slowed by tidal forces. In either 
case, the result would be stronger diurnal cycles than 
those considered in this paper. 

Finally, many directly-imaged planets may have zero 
obliquity and hence no obliquity seasons. We have so far 
focused on true Earth analogs, terrestrial planets orbit- 
ing in the HZ of Sun-like stars. But it is worth noting 
that, due the compact inner- working angle of an interfer- 
ometric direct-imaging mission, many of its target stars 
will be K- and M-dwarfs. Terrestrial planets in the HZ 
of such stars are likely to h ave had their obliquity tidally 
eroded (jHeller et al.ll20li1 ) and therefore exhibit thermal 
phase variations driven only by eccentricity seasons and 
the diurnal cycle. Such lightcurves would be much eas- 
ier to interpret and might allow for a relatively robust 
determination of planetary heat capacity. 

6.3. Constraining the Climate of HZ Exoplanets 

We stated in the introduction that a planet's thermal 
inertia is an attractive "second parameter," after albedo. 
We now expand on this idea, describing how various di- 
rect imaging observations constrain different key climate 
parameters, and how thermal phase measurements fit in 
the big picture. 

0- Dimensional Model: The simplest possible climate 
model balances absorbed stellar flux with thermal emis- 
sion, yielding a planet's emitting temperature. This re- 
quires knowing a planet's albedo, which could be ex- 
tracted from time-averaged thermal observations, but 
not from reflected-light observations because of the 
albedo-radius degeneracy. This is of course insufficient 
to constrain surface conditions and hence habitability. 

1- Dimensional Model: To estimate surface temper- 
ature, one needs to know the pressure of the emit- 
ting layer and surface, as well as the atmospheric lapse 
rate, the rate at which temperature increases as one 
descends in the troposphere. A thermal mission will 
detect any greenhouse gases since by definition they 
absorb in the thermal IR, and could measure sur- 
face pressure via pressure- broadened absorption features 
([Meadows fe Seagerll2010h . An optical mission may be 
able to detect some greenhouse gases in the near-infrared. 
Significantly, a thermal mission can directly probe the 
surface temperature in opacity windows, providing an 
estimate of lapse rate and a check on climate models. 
Optical data will contain, at best, indirect information 
about the planet's temperature. 

Higher- Order Model: To construct a more comprehen- 
sive climate model, many additional planetary charac- 
teristics must be kno wn. An optical m i ssion could con- 
strain rotation rate dPalle et al.l [20081: lOaklev fe Cash! 
20091) and obliquity (Kawah ara fe Fuiiil I2010L 120111 : 
Fujii fe Kawaharall2012| ). Furthermore, the presence of 
oceans could be inferred via their colors (iFord et al.1 
[20011 ICowan et al.ll2009L [2011 iFuiii et al.l 120101b polar- 
izatio n dZugger et al. 1120101 1201 ID and/or specular reflec- 
tion (Willia ms fe Gaidosl I2008i ~ but confounding effects 
include clouds, Robinson et al. 2010, and snow, Cowan 
et al. 2012). These quantities would improve our under- 
standing of planet formation and geophysics, and would 
make it possible to convert observed thermal phase varia- 
tions into estimates of intrinsic seasonal cycles. For plan- 
ets theorized to have zero obliquity, thermal phase vari- 
ations could constrain planetary thermal inertia. But 



meridional heat transport and surface inhomogeneities 
would impede efforts to estimate planetary thermal iner- 
tia on planets with non-zero obliquity. 

An exception for which the planetary thermal inertia 
can be estimated without resorting to periodic climate 
forcing is a planet where the thermal inertia is domi- 
nated by its atmosphere. Recall that the lapse rate is 
a function of surface gravity and specific heat capacity: 
T = — g/c' p , where the effective specific heat capacity of 
the planet, c' p , includes the effects of latent heat. The 
thermal inertia coefficient for a well-mixed atmosphere 
is /i = c' p p s /g, where p s is the surface pressure, and the 
latent heat effects of water vapor are again included im- 
plicitly. A thermal mission of sufficient spectral resolu- 
tion to measure pressure broadening could therefore es- 
timate an exoplanet's atmospheric thermal inertia, pro- 
vided that the planetary mass, and hence surface gravity, 
were otherwise constrained. 

7. CONCLUSIONS 

Each of the three periodic climate forcings — diurnal 
cycle, as well as obliquity and global seasons — contribute 
to the thermal phase variations of habitable zone plan- 
ets. The damping for each of these oscillations is dif- 
ferent. The result is that while time-averaged thermal 
observations will powerfully constrain the climate of di- 
rectly imaged exoplanets, phase variation measurements 
will be very challenging to interpret, unless planets can 
be assured to have zero obliquity. 

In a few limited cases, the lack of thermal variability on 
orbital timescales would provide circumstantial evidence 
of clement surface conditions. On the other hand, very 
large swings in apparent temperature might discount a 
planet's habitability. 

Our major conclusions are: 

1) Meridional heat transport can significantly affect 
obliquity seasons, and varies in relative strength from 
one planet to the next. Surface inhomogeneities in ther- 
mal inertia, meanwhile, lead to a mixing of high ampli- 
tude, small phase-lag signals with low-amplitude, large 
phase-lag signals. These complications make it impos- 
sible to quantitatively estimate a planet's heat capacity 
from obliquity-driven thermal phase variations, even if 
its obliquity is known through other means. 

2) Although the rotation rate of an Earth analog is 
shorter than the expected integration times of foresee- 
able thermal missions, the diurnal cycle shows up in the 
thermal phase variations on orbital timescales, and is the 
dominant effect for tropical observers. For rapidly rotat- 
ing planets, the diurnal cycle and global seasons are both 
damped by heat storage alone, but planetary thermal in- 
ertia varies as a function of forcing frequency, whether 
heat storage operates as a mixed-layer or diffusively. 

For a realistic habitable zone planet with an inhomo- 
geneous surface, non-zero obliquity, and non-zero eccen- 
tricity, thermal phase variations are insufficient for mea- 
suring planetary thermal inertia. 
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